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Abstract. We present a numerical investigation of the cellular burning regime in Type la supernova explosions. This regime 
holds at small scales (i.e. below the Gibson scale), which are unresolved in large-scale Type la supernova simulations. The 
fundamental effects that dominate the flame evolution here are the Landau-Darrieus instability and its nonlinear stabilization, 
leading to a stabilization of the flame in a cellular shape. The flame propagation into quiescent fuel is investigated addressing 
the dependence of the simulation results on the specific parameters of the numerical setup. Furthermore, we investigate the 
flame stability at a range of fuel densities. This is directly connected to the questions of active turbulent combustion (a mech- 
anism of flame destabilization and subsequent self-turbulization) and a deflagration-to-detonation transition of the flame. In 
our simulations we find no substantial destabilization of the flame when propagating into quiescent fuels of densities down to 
~ 10 7 gem -3 , corroborating fundamental assumptions of large-scale SN la explosion models. For these models, however, we 
suggest an increased lower cutoff for the flame propagation velocity to take the cellular burning regime into account. 

Key words. Supernovae: general - Hydrodynamics - Instabilities 

1 . Introduction to the flame thickness it covers more than 1 1 orders of magni- 
tude. Therefore, large-scale simulations provide an insight into 

Type la supernovae (SNe la) are usually attributed to ther- the flame dynamics on the largest scales but have to rely on as- 

monuclear explosions of white dwarf (WD) stars consisting sumptions about me physics of flame propagation on smaller 

03 ■ of carbon and oxygen. Throughout this paper, we will refer scales? which is poorly underst ood so far. However, the mo- 

to the currently favored model in which the combustion pro- tivation for this approach is that the key feature of the SN la 

ceeds in form of a flame starting near the center of the WD and explosion is turbulent combustion. Turbulence is driven from 

traveling outward. More specifically, the initial mode of flame large . SC ale instabilities— mainly the buoyant Rayleigh-Taylor 

propagation is assumed to be the so-called deflagration mecha- instability and the Kelvin-Helmholtz (shear) instability— that 

nism, in which the combustion wave is mediated by microphys- give rfse tQ me formation of a turbu lent eddy cascade. The tur- 

ical transport processes giving rise to a subsonic flame speed. bulent eddies wrinkle the flame and increase its surface . This 

Underlying our assumptions is the canonical single-degenerate enhances the net burning rate and accelerates the flame. In this 

Chandrasekhar-mass SN la model. For a review of SNJa ex- way the energy production rate C an reach values that are suf- 

plosion models we refer to| Hillebrandt & Niemeyer| (|2000j). fident tQ pQwer a SN Ia eX pi osion . The underlying assump- 

Despite the success of recent a ttempts to numerically tion of this model is that the flame evolution is dominated 

model Type la super nova explosions jReinecke et all l2002bt by the turbulent cascade originating from instabilities on large 

iGamezo etafl l2003h . some basic questions regarding the ex- scales. Although there exist some theoretical ideas corroborat- 

plosion mechanism remain unanswered. This is mainly due to ln g the assumption of flame stability on small scales, a justifi- 

the fact that the whole range of scales relevant to this problem cation by means of a hydrody namical simulation has not pro- 

cannot be resolved in a single simulation in the foreseeable fu- vided definitive answers ye t dNiemever & Hillebrandl Il995t 

ture. From the radius of the exploding white dwarf star down iNiemever & WooslevL| 19971 The results that will be reported 

in the present paper are, however, a contribution into this direc- 

Send offprint requests to: F. K. Ropke tion. 
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The work presented in the following stands in direct succes- 
sion o f Rop ke et all fe003l) . where the basic method was intro- 
duced and applied to describe flame propagation into quiescent 
fuel for a density of the unburnt material of p„ = 5xl0 7 g crrT 3 . 
This example demonstrated that the applied method appropri- 
ately models the flame propagation on specific scales in the SN 
la explosion. Moreover, it confirmed for the first time by means 
of a full hydrodynamical simulation that the stabilization of the 
flame front in a cellular pattern known from chemical flames 
holds for thermonuclear flames under the conditions of SN la 
explosions. In the present paper we will extend the study of 
flame propagation to a wider range of fuel densities. The aim is 
to test the stability of the cellular pattern in dependence of this 
parameter. This will provide an overview over possible effects 
of flame propagation into quiescent fuel. We will additionally 
investigate the flow field resulting from the flame configuration 
in greater detail and measure the effective flame propagation 
velocities for the resulting flame structures. 

The theoretical context of this study will be introduced in 
Sect. 13 After a brief description of the applied numerical meth- 
ods in Sect. [3] we will present the results of numerical simula- 
tions of the flame evolution in Sects. @] and [5] Finally, conclu- 
sions regarding the significance of the numerical investigations 
for large-scale SN la models will be drawn. 

2. Theoretical background 

Our study of flame evolution focuses on scales where flame 
propagation decouples from the turbulent cascade. On larger 
scales, this cascade dominates flame propagation by wrinkling 
the flame and enhancing its surface. But in order to contribute 
to the flame wrinkling, a turbulent eddy has to deform the flame 
significantly in a time shorter than the flame crossing time over 
the spatial extent of that eddy. Thus, equating the eddy turnover 
time to this time scale, one obtains the cutoff scale /cibs 



v'fcbs) = si, 



(1) 



where si denotes the so-called laminar burning velocity of a de- 
flagration flame, i.e. the propagation speed of a planar deflagra- 
tion wave which is determined by microphysical transport. The 
turbulent velocity fluctuations at length scale / are represented 
by v'(0- Equation Q defines a length scale labs, whi ch was 
identified with the Gibson scale for chemical flames (Peters, 
119861) . We will make use of this term in the context of ther- 
monuclear flames in SN la explosions, too. 

In our model, we describe the flame as a discontinuity be- 
tween burnt and unburnt states ignoring any internal structure. 
This simplification is exploited by the numerical method of 
flame tracking applied in our simulations, resulting in an in- 
creased dynamical range compared to models that fully resolve 
the flame structure. In order to justify our thin flame approxi- 
mation, we provide a rough estimate of the Gibson scale. The 
turbulent velocity fluctuations V in eq. Q are given by the 
scaling law of the turbulent eddy cascade. Kolmogorov scal- 
ing, v'(l) oc Z 1 / 3 , provides a reasonable approximation. Hence 
one obtains from 

,3 



fcibs - L 



v'(L) 



The large scale instabilities give rise to the formation of as- 
cending bubbles. A typical size of a Rayleigh-Taylor bubble is 
L ~ 10 7 cm, which defines the integral scale of the turbulence. 
The rise of these bubbles invokes shear velocities v'(L) of about 
10 7 cms _1 . Values for the laminar burning velocity have been 
determined by means of one-d imensional simulations resolving 
the flame structure bv lTimmes & Wooslevl i ll 9921) . The authors 
also provide the following fitting formula to the data which we 
will use in our simulations: 



si = 92.0 x 10 : 



\2x 10 9 / 



X( lz C) 



0.5 



0.889 



cms 



(3) 



However, we note that the fit does not reproduce the simulation 
results as accurately as the authors claim for the fuel density 
range that will be investigated in the following. According to 
this formula, the laminar burning velocity at a fuel density of 
5 x 1 7 g cm -3 is of the order of 1 6 cm s~ 1 . This yields a Gibson 
scale of ~10 4 cm. Thus, keeping in mind the indicative charac- 
ter of our study, we will arrange our simulations at length scales 
around this value, which is still well-separated from the flame 
width. 

Below the Gibson scale, flame propagagation does not in- 
teract with the turbulent cascade, a s in the context of SN 
la exp losion was first pointed out by iNiemever & Hillebrandtl 
(1995). Here, it would propagate as a laminar flame. However, 
it is well-known that laminar flames are subject to a hydrody- 
namical instability (Darrieusl Il938l Landau! Il944 . The ori- 
gin of this so-called Landau-Darrieus instability (LD instabil- 
ity henceforth) is the refraction of the streamlines of the flow 
on the density contrast across the flame. Mass flux conservation 
leads to a broadening of the flow tubes in the vicinity of bulges 
of the perturbation. Here the local fluid velocity decreases and 
becomes lower than the fluid velocities at +oo. These velocities 
correspond to the laminar burning speeds with respect to fuel 
and ashes in the frame of reference comoving with the flame. 
Therefore the burning velocity s/ of the flame is higher than 
the corresponding local fluid velocity and this leads to an in- 
crement of the bulge. The opposite holds for recesses of the 
perturbed front. In this way the perturbation keeps growing. 

This can be observed in our simulations. Figure^ depicts 
the absolute value of the fluid velocity normalized to the lami- 
nar burning velocity of the flame in the linear regime of flame 
evolution in one of our simulations. It is well in agreement with 
the theoretical expectations and confirmes that our numerical 
method is capable to reproduce the LD instability. This is not 
trivial as can be seen from the failure of th e passive implemen - 
tation to develop the correct flow field (cf. Rop ke et al ]l2003h . 



By means of a linear stability analysis. lLandaul \ 19441) de- 
rived a dispersion relation between the growth rate of the am- 
plitude a> and the wavenumber k of the perturbation: 



wld = ksi 



1 + p 



(4) 



(2) 



with/7 =p u /Pb- 

This instability, however, is not necessarily a contradiction 
to the assumption of flame stability on small scales. In fact, 
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(a) 




(b) 

Fig. 1. Fluid velocity in the simulation with p u = 5 x 10 7 g cm" 3 
and 100 x 100 cells resolution. The x-axis and the y-axis show 
position in cm; the z-axis represents \v\/s/ (a) snapshot at t — 
0.01 s, (b) snapshot at t = 0.09 s 



one observes stable flame propagation in chemical combus- 
tion experiments. The reason for this is that the flame stabi- 
lizes in a cellular pattern as soon as the growth of the perturba- 
tion enter s the nonli near regime. This effect has been explained 
bv lZerdovichl ll966) in a geometrical approach. The propaga- 
tion of an initially sinusoidally perturbed flame is followed by 
mean s of Huygens' principle (see also Fig. 1 in R opke etall 
120031) . At the intersections of two bulges of the front a cusp 
forms after a while. Here Huygens' principle is no longer ap- 
plicable and the flame evolution enters the nonlinear regime. 
The propagation velocity at the cusp exceeds sc. 



where 9 denotes the inclination angle of the front adjacent to 
the cusp. This effect balances the perturbation growth due to 
the LD instability and leads to a stabilization of the flame in 
a cellular configuration. Figure^) shows the absolute value of 
the velocity after a cusp has formed. 

Assuming a parabolic shape of the cells, the temporal evo- 
lutio n of the am plitude of the perturbation x can now be written 
as dZerdoviclill966l) 



— = U) LD X -k SIX 

at n z 



(6) 



The quadratic damping term accounting for the cellular stabi- 
lization introduces an extension to the linear stability analysis 
resulting in the dispersion relation (j4ji. If one further assumes 
that the flame front advances with the constant laminar burning 
speed si independently of the flame shape, the effective speed 
v e tf of the mean position of a flame wrinkled due to instabilities 
is given in terms of its surface A' compared to the surface of a 
planar flame A(0): 



A' 



Veff = 



(7) 



The flame thus accelerates with increasing surface area 
and the task to determine the effective flame propagation ve- 
locity reduces to the determination of the increase in flame 
surface area. From the assumption of a parabolic cell shape 
IZel'dovich et alJ ([I58O) derived the acceleration of the flame 
corresponding to the increase in its surface: 



Vcell . 

= 1 + e, 

Si 



e= n -ll- 1 - 

24 



(8) 



si 
cos 8' 



(5) 



The cellular stabilization is the reason that large-scale 
SN la simulations assume flame stability on small scales. 
This is, however, not yet well confirmed. The cellular pat- 
tern may break up under certain conditions, e.g. density 
jum ps over the flame fron t or interaction with turbulence. 
iBlinnikov & Sasorovl dl996h suspected that a fractalization of 
the cellular flame may drastically accelerat e the flame at low 
fuel densities and iNiemever & Hillebra ndt ( 1995) reported on 
indications that the flame may destabilize at fuel densities 
of p u a; 5 x 10 7 gcrrT 3 . However, a failure of the cellu- 
lar flame stabilizat ion at p , . = 5 x 10 7 gcnT 3 could not be 
confirmed by Ropke et al. (2003). Nevertheless, the possible 
loss of flame stability at certain densities of unburnt material 
could have drastic impact on large-scale SN la models and 
has to be soundly tested. Besides the break-down of one fun- 
damental assumption of these models, it could also have ef- 
fects that seem appealing. Spectra produced by empirical one- 
dimensional SN la models fit the observations remarkably well 
if a transition of the flame propagation mode from subsonic 
deflagrat ion to supersonic detona t ion is applied at low den- 
sities faoflich & Khokhlovt Il996l llwamoto et all Il999h . The 
drawback of these delayed detonation mode ls is that th e mech - 
anism providing the transition is unclear dNiemeverl [l999). 
As one possibility, an active production of turbulence and 
self-turbulization of the flame after a break-down of its cel- 
lular stabilization was proposed under the name active turbu- 
lent combustion by Niemever & WoosleyJ (119971) based on an 
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idea by Kerstein ( 1996). Following the reasoning provided by 
lNiemever"&W)oslevl ( 119971) . this effect is directly connected 
to the burning regime we are aiming at in our simulations. 

The above considerations point out the need for a thorough 
investigation of the cellular burning regime in SN la explosions 
at scales around or below the Gibson length. Numerical simu- 
lations addressing this question will be presented in Sects. |4] 
and|5] 

3. Numerical methods and simulation setup 

The numerical methods we apply are based on the work by 
iReinecke et al.1 d 19991) and have been described by R opke et all 
l l2003l) . Therefore we will be brief in this and on ly mention 
some keywords here refering to Rop ke etafl J200 3) for the de- 
tails. 

Since our simulations aim on the range around the Gibson 
scale, which is well above the flame thickness for the fuel den- 
sities we are going to apply, it is justified to model the flame 
as a discontinuity between unburnt and burnt states. This flame 
description ignores any internal structure and therefore does 
not intrinsically provide the value of laminar burning velocity. 
This also implies that, in a first approach, hydrodynamics and 
flame propagation can be modeled in separate steps. 

The hydrodynamics is modeled applying the piecewise 
parabolic method (PPM) — a highe r order Godunov scheme de- 
veloped by Colella & Woodward ( 1984) — in the Prometheus 
implementation (iFrvxell et all Il989l) . The corresponding part 
of the code treats only the fluid dynamics and does not ac- 
count for the flame propagation. Since the flame is modeled 
in the discontinuity approximation, its propagat ion can be de - 
scribed applying the level-set method dOsher&Sethiani ri988). 
Rop ke et al.1 (12003) could show that the key feature that enables 
us to reproduce the theoretically anticipated flame evolution 
resulting from the LD instability is a specific implementation 
of flame/flow coupling, namely t he in-cell reconstru c tion/fl ux- 
splitting technique proposed by Smil ianovski et alJ dl997l) . It 
couples the flame propagation accurately to the flow and there- 
fore allows the investigation of effects that originate from hy- 
drodynamics, such as the LD instability. Our implementation 
corresponds to what was termed "complete imp lementation of 
the level-set method" by Reinec ke et alJ ([1999) with some mi- 
nor changes (cf. lRopke et alJl2003l) . All simulations presented 
in the following were performed in two spatial dimensions. 

The basic approach is to study the flame evolution under the 
influence of a small perturbation (which can be introduced by 
noise or — as in our case — be initially imprinted to the flame 
structure). This perturbation should grow due to the LD in- 
stability. In the nonlinear regime, the formation of a cellular 
pattern should inhibit the further growth of the perturbation 
and s tabilize the flame. In a previous publication (Ropke et al., 
120031) we demonstrated that this mechanism holds for ther- 
monuclear flames in degenerate matter. From theoretical con- 
siderations and semi-analytical models (resulti ng for instance 
from the Sivashinsky equation, e.g. iGutman & Sivashinskvl 
Il990l) . it can be expected that the general features of the evo- 
lution of the flame front shape depend significantly on the nu- 
merical setup used in the simulations. Possible parameters that 



influence the flame evolution are the overall geometry of flame 
propagation, the width of the computational domain compared 
to the length scale of perturbations, the resolution, the bound- 
ary conditions, sources of numerical noise (which is likely to be 
different when the simulation is parallelized to a varying num- 
ber of sub-processes), and the initial flame shape. Therefore 
care has to be taken in choosing the specific setup depending 
on the questions that are addressed by the simulations as well 
as in the conclusions drawn from the simulation results. 

Regarding the overall flame geometry, two cases are com- 
monly studied in the literature: an on average planar flame ge- 
ometry and a flame that is on average circularly expanding. 
Although a naive approach would choose the second scenario 
for the supernova explosion, it is probably not the appropriate 
description of the flame propagation there. The flame evolu- 
tion on scales of the star, where expansion effects are most pro- 
nounced, quickly becomes dominated by the Rayleigh-Taylor 
instability. This leads to a flame evolution completely diverg- 
ing f rom a circular (or in three dimensions spherical) geometry 
(cf. lReinecke et al.l2002al) . It rather proceeds in raising bubbles 
of burnt material. Nevertheless, it could be argued that these 
structures again partly resemble a spherically expanding geom- 
etry. Although the case of a circularly expanding flame reveals 
very interesting physical effects (as repeated mode splitting of 
the cells resulting from expansion effects and a possibly result- 
ing fr actalization of the flame front, cf. Blinni kov & Sasorovl 
1996j), we will not follow this approach here for two reasons: 
First, in the scope of our implementation it is prohibitively ex- 
pensive to follow a circular flame evolution for a sufficiently 
long time and, second, we aim on effects on scales around the 
Gibson length, at which global expansion effects are negligible. 
This has a technical advantage. It is much simpler to keep an 
overall planar flame in the center of the domain by choosing a 
comoving frame of reference. This enables us to study the long 
term flame evolution without requiring large computational do- 
mains. Thus it becomes possible simulate the long-term flame 
evolution, as will be presented below. Note that the choice of an 
overall planar flame geometry determines a distinct flame evo- 
lution and possibly excludes the mechanism of repeated mode 
splitting. 

The influence of the width of the computational domain has 

■ 1 1 ' 1 

already been discussed by Ropke et al. ( 2003). The simulations 
presented there lead to the conclusion that in case of sufficient 
numerical resolution the flame stabilizes in a single domain- 
filling cusp-like structure for periodic boundaries transverse to 
the direction of flame propagation. However, with increasing 
domain width (corresponding to a higher resolution in the sim- 
ulations) this fundamental structure became superimposed by a 
short-wavelength cellular pattern. Thi s is in accord wi th results 
from semi-analytical studies dGutman & Sivashinskvl 1 1990). 

The discussion of the influence of the boundary conditions 
and the initial shape of flame perturbations will be postponed 
to the next section. 

The "experimental setup" applied our simulations is the fol- 
lowing: The spatial extent of the computational domain was 
set to correspond to scales around the Gibson length and was 
fixed according to the discussion above. The flame was initial- 
ized in the center of the domain in an on average planar ver- 
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Fig. 2. Flame propagation into quiescent fuel at p„ — 5 x 10 7 g cirT 3 ; resolution: 200 x 200 cells, periodic boundary conditions; 
Snapshots taken at (a) 0, (b) 5, (c) 7.5, (d) 10, (e) 15, and (f) 30 growth times r LD = a>~^ of a perturbation with A = 200Ax. The 
vorticity £ (cf. eq. l|9}) is color-coded and the flame position is indicated by solid white curves. 
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Fig. 3. Flame propagation into quiescent fuel at p u — 5 x 10 7 g cirT 3 ; resolution: 200 x 200 cells, reflecting boundary conditions; 
Snapshots taken at (a) 0, (b) 5, (c) 7.5, (d) 10, (e) 15, and (f) 30 growth times r LD = a>~^ of a perturbation with A = 200Ax. The 
vorticity £ (cf. eq. l|9}) is color-coded and the flame position is indicated by solid white curves. 
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tical shape with the unburnt material on the right hand side 
and the burnt material on the left hand side, so that in lab- 
oratory frame of reference the flame would propagate to the 
right. The domain was set to be periodic in y-direction. On the 
left boundary of the domain an outflow condition was enforced 
and on the right boundary we imposed an inflow condition with 
the unburnt material entering with the laminar burning veloc- 
ity si. This would yield a computational grid comoving with 
a planar flame. However, the LD instability leads to a growth 
of the perturbation and therefore increases the flame surface. 
According to eq. Q, this causes an acceleration of the flame 
and therefore it is necessary to take additional measures in or- 
der to keep the mean position of the flame centered in the do- 
main. One possibility is to detect the mean location of the flame 
and to simply shift the grid to keep the flame in the center. This 
method is consistent with the boundary conditions applied in 
jc-direction. The described method allows the study of the long 
term flame evolution. In order to induce the development of 
perturbations, we usually imposed a sinusoidal perturbation on 
the initial flame shape. 

The state variables were set up with values for the burnt and 
unburnt states obtained from (pseudo-)one-dimensional sim- 
ulations performed with the "passive implementation" of the 
level-set method ( Reinecke et all ll999). imposing a value for 
the density of the fuel. A compilation of the relevant setup val- 
ues for different fuel densities is given in Table [2 Note that 
frequently in this paper fuel densities refer to a label of a spe- 
cific set of values rather than giving the accurate fuel density. 
The laminar burning speed is calculated according to eq. 0. 

4. General features of flame evolution 

The flame evolution for the exemplar y case of a fuel d ensity 
of 5 x 10 7 gcirr 3 was described by Rop ke et alJ (12003). The 
results obtained in this reference were 

- In the linear regime of flame evolution the initial perturba- 
tion grows due to the LD instability. 

- The formation of cusps at recesses of the flame front leads 
to a stabilization of the flame structure in a cellular shape 
in the nonlinear stage of evolution. 

- Different resolutions were tested. The flame stabilizes in a 
single domain-filling cusp-like structure for resolutions of 
100 x 100 cells and above. 

- In highly resolved simulations, the fundamental cusp-like 
structure became superimposed by a short-wavelength cel- 
lular pattern. This effect is in agreement with semianalytical 
results modeling the flame evolution with the Sivashinsky- 
equation dGutman & SivashinskvHl99oT) . 

- The growth of the perturbation amplitude in the linear 
regime is consistent with eq. |4] for sufficient resolution 
(above [200 cells] 2 to [300 cells] 2 ). 

- The increase in flame surface area agrees reasonably well 
with the gain in the effective flame propagation speed. 

In this section we extend this investigation by presenting 
two simulations of the long-term flame evolution, focusing on 
the nonlinear regime. The two simulations were carried out on 
a regular Cartesian grid with 200 x 200 cells with a cell width 




8 x 10 3 cm 



Fig. 4. Evolution of the flame front for p„ = 5 x 10 7 gcnT 3 
(a) with periodic boundaries in v-direction, (b) with reflect- 
ing boundaries in y-direction. Each contour marks a time step 
of 0.5 TLD(/lpert = 3.2 x 10 4 cm). The contours are artificially 
shifted for better visibility and do not reflect the actual flame 
propagation velocity. 

of Ax = 160cm. The fuel density was 5 x 10 7 gem" 3 . In both 
cases we perturbed the flame front initially in a sinusoidal way 
with a wavelength of 1/6 of the domain width, but we applied 
different boundaries transverse to the flame propagation. 

Figure |3 illustrates the flame evolution in case of periodic 
boundary conditions in y-direction. The first snapshot shows 
the initial setup with the sinusoidal perturbation imposed on the 
flame. In the second snapshot the flame has already reached the 
nonlinear regime and exhibits a cellular shape. Here, the onset 
of an irregularity of the pattern is apparent. This effect leads to a 
merging of cells observable in the next three snapshots. Finally, 
the flame reaches a steady state which is a single domain-filling 
cusp. Note that the snapshots in Fig. |2 were not taken at equal 
time intervals. The final flame configuration was followed for 
more than 10 growth times tu> = u> L l D corresponding to a per- 
turbation of the width of the domain in order to ensure its sta- 
bility (cf. Fig. Eli). 

In the snapshots of the flame evolution depicted in Fig. |2] 
the vorticity of the flow field, £ = V x v, is color-coded. For our 
two-dimensional simulations the (scalar) vorticity reads 



Note that we color-code the logarithm of the absolute value of 
the vorticity and therefore even slight deviations from zero vor- 
ticity are visible. The flame produces vorticity in the burnt ma- 
terial. This effect is especially strong behind cusps. Analytical 
approximations of the flame evolution frequently make the as- 
sumption of a potential flow ahead of the flame. In contrast to 
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Table 1. Setup values for the simulations of the flame evolution. 



Label 


p„[gcm 3 ] 


p 6 [gcm 3 ] 




At 


Jilcms" 1 ] (TW) 


e,.„[ergg '] 


ey,[ergg '] 


1 x 10 7 


9.969 x 10 6 


2.691 x 10 s 


3.70 


0.575 


2.39 x 10 5 


1.619 x 10 17 


8.59 x 10 17 


1.25 x 10 7 


1.232 x 10 7 


3.830 x 10 6 


3.22 


0.526 


2.86 x 10 5 


2.005 x 10 17 


8.99 x 10 17 


2.5 x 10 7 


2.484 x 10 7 


8.69 x 10 6 


2.86 


0.482 


5.01 x 10 5 


2.56 x 10 17 


9.50 x 10 17 


5x 10 7 


4.988 x 10 7 


2.071 x 10 7 


2.41 


0.413 


8.74 x 10 5 


3.584 x 10 17 


1.051 x 10 18 


7.5 x 10 7 


7.50 x 10 7 


3.345 x 10 7 


2.24 


0.383 


1.21 x 10 6 


4.29 x 10 17 


1.13 x 10 18 


1 x 10 8 


1.00 x 10 8 


4.789 x 10 7 


2.09 


0.352 


1.53 x 10 s 


4.89 x 10 17 


1.182 x 10 18 


1 x 10 9 


1.00 x 10" 


6.80 x 10 8 


1.47 


0.190 


9.75 x 10 6 


1.27 x 10 18 


1.96 x 10 18 



that, we observe the development of vorticity in the fuel re- 
gion. This, however, is not surprising, since the flame propa- 
gates subsonically and the regions downstream and upstream 
of the flame are causally connected. 

The merging proceeds in a way that some of the cells grow 
and the smaller cells disappear in the cusp between larger cells 
as can be clearly seen in Fig. |4t- This effect forms the ba - 
sis of flame stability as explained by Zerdovic h et al 1 ill 98(1 
There, the authors addressed the phenomenon that flames in 
experiments are often observed to stabilize in large cells. This, 
however, seems to be in contradiction with the result from 
Landau's linear stability analysis (eq. @), because it predicts 
higher growth rates for perturbations with smaller wavelengths, 
which thus should dominate the flame structure. Even taking 
into account a modified dispers ion rel ation resulting from a fi- 
nite flame thickness jMarksteinlll95ll) cannot explain this phe- 
nomenon. |ZeFdovicl^^nTI^I3rarg ue that the stable long- 
wavelength flame pattern is a result of the flow that establishes 
upstream of the flame. A velocity component tangential to the 
flame is directed toward the cusp and advects small perturba- 
tions into this direction. A gradient in this tangential velocity 
stretches the perturbation wavelength and thereby retards its 
growth. Fi nally, the small p erturbation disappears in the cusp 
('originally IZel'dovich et al.l Jl980l) analyzed the flame propa- 
gation in tubes, but they point o ut the similarity of this con- 
figuration with a cellular flamel. IZel'dovich et alJ (| l98Q!) give 
an analytical description of the effect applying a WKB-like ap- 
proximation. 

The described mechanism acts also in our numerical model, 
as was discussed by Ro oke et al.l (120031) in connection with the 
flame structure superposed by a short-wavelength cellular pat- 
tern for high-resolved simulations. A close-up of the flow field 
in vicinity of a cusp taken from one of our simulations (see 
Fig. |5} demonstrates the ability of our implementation to re- 
produce the effect. The arrows represent the velocity field after 
the formation of a cusp. Note that the streamlines converge to- 
ward the cusp ahead of the front. This leads to the formation of 
a layer upstream of the front in which the fluid velocity is di- 
rected toward the cusp. Thus the required flow field providing 
stability of long-wavelength patterns establishes in our simu- 
lations. In the downstream region, the flow diverges from the 
cusp. 




o l.o xlo 4 2.0 x to 4 3.0 xlo 4 

t x [cm] 

Fig. 5. Velocity field in a simulation with a fuel density of p u = 
5 x 10 7 g cirT 3 and 100 x 100 cells resolution at t = 0.09 s; solid 
curve indicates flame position. 

The flame evolution in case of reflecting boundary condi- 
tions in y-direction is shown in Fig. [5] The mechanism of cellu- 
lar stabilization acts similar to the case of periodic boundaries. 
However, the alignment of the steady-state of the flame shape 
is different. Whereas the periodic case developed a cusp in the 
center of the domain, in the reflecting case the crest of the pat- 
tern centers in the domain and the recesses of the front form a 
structure similar to a "half-cusp" at the boun daries. This case 
is ana logous to the configuration studied by IZel'dovich et alJ 
(1980). 

Another aspect of these two simulations should be noted. 
Ropke et al. 1 d2003l) studied the evolution of a flame that was 
initially perturbed by only one single domain-filling mode. 
This was done in order to be more efficient in the investiga- 
tion of the linear growth due to the LD instability and it was 
argued that a single domain-filling structure establishes inde- 
pendently of the exact shape of the initial perturbation. This 
is reconfirmed by the presented simulations in accord with 
expectations from numerical simulations of the Sivashinsky- 
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3.2 X LO 4 





S x 10 J tan 



S x LO 3 ran 




8.0 x LO 3 




8 x LO 3 cm 



S x LO 3 ran 



-3 



Fig. 6. Evolution of the flame front at (a) p u — 2.5 x 10 7 g cm 4 Fig. 7. Evolution of the flame front at (a) p„ = 1 x 10 7 g cm 
and (b) p„ = 5 x 10 7 g cm 4 . The contours mark evolution steps and (b) p u = 1.25 x 10 7 g cirT 3 . The contours mark evolution 
of 0.2 tld- Every 10th conture is highlighted bold for better steps of 0.2 tld- Every 10th conture is highlighted bold for 



visibility. 



better visibility. 



equation (e.g. lGutman & Sivashinskvll99C ) and from the pole - 
decomposition solutions of that equation (Thu al et all 11985). 
In both the simulation with periodic boundaries (cf. Fig. [2) and 
the simulation with reflecting boundary conditions (see Fig. [3 
the cells merge. The steady state of the flame resulting from 
this process is a single domain-filling cell. 



5. Flame stability at different fuel densities 

Further numerical experiments addressed the flame stability at 
different fuel densities. With lower fuel density, the Atwood 
number — defined as At = \p u -pb\/(Pu+Pb) — and, equivalently, 
the density contrast p. over the flame front increase (cf. Table 
At the same time the laminar burning velocity of the flame 
decreases. In the following, tests of the flame stability at a va- 
riety of fuel densities will be presented. These use the setup 
introduced in Sect. |3]with periodic boundaries in the direction 
transverse to the flame propagation. We again refer to Tabled 
for a compilation of the setup values. The resolution chosen for 
this study was 200 x 200 cells. From the simulations presented 
in the previous section, it follows that an initial flame perturba- 
tion with a wavelength corresponding to the width of the com- 
putational domain is sufficient in order to see the flame stabi- 
lize in a universal steady state. In the case under consideration 
one expects a single centered domain-filling cusp-like struc- 
ture. Although the restriction to an initial domain-filling per- 
turbation wavelength may appear rather artificial, it avoids the 
time-consuming cell merging process before the final steady- 
state flame structure is reached. 



Figure[6]shows the temporal evolution of the flame front for 
fuel densities of 2.5 x 10 7 g cnT 3 and 5 x 10 7 g cm 4 . In both 
cases the flame stabilizes in a single domain-filling cusp-like 
structure. Note, that the evolution time has been normalized to 
the growth time of the LD instability tld = corresponding 
to the initial perturbation wavenumber. This ensures compara- 
bility between the different simulations. As in Fig.|4] the spac- 
ing between individual contours is, however, artificial since the 
original simulations were performed in the frame of reference 
comoving with the flame. It has been chosen in a way that the 
contours fill the plot window and does not reflect flame propa- 
gation. 

The outcome from simulations with even lower fuel den- 
sities reveals a completely different flame evolution. Figure 
shows the flame for fuel densities of 1 x 10 7 g cm 4 and 1 .25 x 
10 7 g cm 4 . Although the initial flame evolution resembles that 
of p„ « 5 x 10 7 g cm 4 , the forming cusp lacks long-term sta- 
bility. As can be seen from the plots, the stable structure breaks 
up from the cusps outward and the flame subsequently evolves 
in an irregular pattern. The origin of these effects is most likely 
numerical. For a fuel density of 1 x 10 7 gcm 4 a simulation 
run with 600 x 600 cells resolution did not result in a disrup- 
tion of the cusp. Thus, a high grid resolution is required in 
order to describe a stable cusp properly at low fuel densities 
(< 1.25 x 10 7 gem 4 ). 

A rather unexpected result is obtained from simulations for 
fuel densities around 1 x 10 8 g cm 4 . Here, the flame does not 
propagate in a stable manner in the first stages of the flame 
evolution. This is illustrated by Fig. [H] Small-wavelength per- 
turbations superpose the initial flame shape shortly after the be- 
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ginning of the simulation. These appear to dominate the flame 
evolution for a while (see in particular Fig. |8j)) but then the 
flame stabilizes in the single-cell structure. This implies that 
(i) the flame structure is less stable against small-wavelength 
perturbations at these fuel densities and (ii) the mechanism for 
a stabilization in a preferred long-wavelength pattern as ana- 
lytically predicted (see the discussion in Sect. |4} does finally 
stabilize the flame evolution. 



IX LO 4 



Table 2. Flame propagation velocities in the cellular regime. 




S x LO 3 cm 




8 x LO 3 cm 



Fig. 8. Evolution of the flame front at (a) p u — 7.5 x 10 7 g cm' 3 
and (b) p„ = 1 x 10 8 g cirT 3 . The contours mark evolution steps 
of 0.2 Tld- Every 10th conture is highlighted bold for better 
visibility. 



The initial flame destabilization becomes less pronounced 
in flame evolution at even higher densities. Figure |9] shows an 
example with p u — 1 x 10 y g cirT 3 . 



Pu [gem 3 ] 




Vcell/ Si 




(measured) 


(according to eq. {S}) 


2.5 xlO 7 


1.23 


1.17 


5 xlO 7 


1.20 


1.14 


7.5 xlO 7 


1.25 


1.13 


1 xlO 8 


1.30 


1.11 


1.25 xlO 8 


1.35 


1.10 


1 xlO" 


1.09 


1.04 



We measured the flame acceleration in these simulations 
via the increase in flam e surface area (cf. e q. 0). This method 
has been discussed by Roo ke et alJ (|2(}03). Although there are 
some minor differences compared to measuring the effective 
flame propagation velocity directly (most important, eq. @ 
strictly holds only in case of a flame idealized as a discontinu- 
ity, which is fulfilled to a reasonable degre e in our sim ulations, 
but probably not quite exactly, cf. iRopke et alJl2003() . this ap- 
proach can be expected to yield rather reliable results. Some 
examples of the temporal evolution of the flame surface area 
are plotted in Fig.^| This figure illustrates the different cases 
of flame evolution. While the simulations of low fuel densi- 
ties (Fig. E5|i) fail to converge, the flame reaches a steady-state 
propagation velocity for higher fuel densities. The initial desta- 
bilization of the flame pattern and final stabilization for fuel 
densities around 10 8 g cirT 3 is apparent in the evolution of the 
corresponding flame surface area (Fig.llOt). 

Table [2] compares the final steady-state velocities of the 
cellular flames to the theoretical predictions from eqs. l|8}. 
These equations, of course, ignore any structure superimposed 
to the front the flame front and additionally assumes a parabolic 
shape of the smooth cells. Both effects obscure the expected 
trent of decreasing v ce u/s/ with increasing fuel density. This can 
be clearly seen from Figs. |6j]9] where the deviations are most 
pronounced for fuel densities around 10 8 g cm 3 . Therefore the 
measurements agree only in the order of magnitude with the 
theoretical values. This indicates that ignoring the superim- 
posed small-wavelength pattern is a too simple approach and 
underestimates the actual flame surface. 




8 x L0 3 cm 



Fig. 9. Evolution of the flame front at p„ — 1 x 10 9 g cirT 3 . The 
contours mark evolution steps of 0.2 r LD . Every 10th conture 
is highlighted bold for better visibility. 



6. Conclusions 

The simulations presented in Sects. |4] and [5]provide a deeper 
insight into the dynamics of thermonuclear flame fronts under 
conditions of SN la explosions. One fundamental assumption 
was that the flame propagates into quiescent fuel. 

In Sect. |2 we addressed the question of the influence of the 
boundary conditions and the shape of the initial flame pertur- 
bation on the long-term flame evolution. The conclusions that 
can be drawn from our simulations are: 

1 . In the given setup, the initial perturbation shape is not re- 
tained. Although at first it develops into a cellular struc- 
ture of the same wavelength in the nonlinear regime, the 
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(a) 



Fig. 10. Increase in flame surface area at fuel densities of (a) p u 
1 x 10 s g cm" 3 , (d) Pu = 1 x 10 9 g cm" 3 . 




(b) 



(d) 

1.25 x 10 7 gcm- 3 , (b) p u = 5 x 10 7 gcirT 3 , (c) p u 



flame later shows the tendency to align in the shape of 
a single domain-filling cell. The transition to this steady- 
state pattern happens by growth of preferred cells from 
the initial perturbation and the disappearance of smaller 
cells in newly formed cusps. This effect gives rise to a 
"mer ging process" of cells, consisten t with numerical solu- 
tions (Gutman & Sivashinskv, 1990) and with pole decom- 
position analysis ijT hual et alj 119851) of the Sivashinsky- 
equation ( Sivashinskv! 1 1 977h 

2. The merging of the small cells proceeds smoothly without 
drastic effects on the flame shape. In particular it does not 
lead to a loss of flame stability. 

3. The cases of reflecting and periodic boundary conditions 
differ in the alignment of the ultimately developed cell. In 
the case of periodic boundary conditions the cusp devel- 
ops at the center of the domain and the crest tends toward 
the boundary. In principle no specific alignment should be 
preferred for periodic boundaries. This is not true for re- 
flecting boundary conditions. Here, the crest centers in the 
domain, simila r to the simulation with adiabatic boundaries 
presented by Gutma n~& Sivashinskv! dl990h . 

Thes e results support the generality of the simulations pre- 
sented by Rop ke et al and were used to set-up the study 
of the influence of the fuel density on the flame stability that 



was presented in Sect. [5] He re, the general featu res of the flame 
evolution resemble those o f iRopke et alJ §003) and Sect.E|of 
the present paper. However, some peculiarities were observed. 
At very low fuel densities (< 1.25 x 10 7 gem -3 ), the cusp-like 
flame shape breaks up. It was argued that this is a numerical ef- 
fect. Nevertheless, if the reason for the disruption of the stable 
flame pattern is numerical noise, then one could conclude that 
flames at lower fuel densities are likely to become more sensi- 
tive to velocity fluctuations in the fuel region. These, however, 
can be expected to be present in realistic scenarios of SN la 
explosions. Therefore we will address that question in a subse- 
quent study. 

A destabilization of the flame structure shortly after initial- 
ization was observed at densities around 10 8 g cm- 3 . However, 
after ~2tld the flame stabilized again forming the cusp-like 
steady-state structure. Thus, we conclude that the flame tends 
to stabilize in a cellular pattern in the parameter space explored 
in our simulations. 

The implications of the presented investigation for large- 
scale SN la models can be summarized as follows: 

(i) The flame is expected to propagate in a stable way 
at least down to fuel densities of around 10 7 gcrrr 3 . This 
corroborates fundamental assumptions of large-scale SN la 
models (e.g. iReinecke et al.ll2002bl) which rely on a smooth 
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flame evolution on unresolved scales. We dis agree with pre- 
vious simulations by Niemever & Hillebrandt ( 1995), who re- 
ported on a possible loss of flame stability at a fuel density of 
5 x 10 7 gcirT 3 . Therefore, a self-turbulization and subsequent 
drastic acceleration of the flame seems unlikely — at least in the 
case of propagation into quiescent fue l. No convincing indi- 
cations for activ e turbu lent combustion dNiemever & Wooslevt 
119971: iKersteiri Il996l) 1 were observed in our simulations. 
However, flame propagation into a turbulent flow field may re- 
veal completely different flame structures, since one expects 
strong enough turbulent eddies to disrupt the cellular flame pat- 
tern. The arguments above provide strong hints into this direc- 
tion. 

(ii) Large-scale S N la simulations as performed by 
iReinecke et alJ J2002b) apply a subgrid-scale model to deter- 
mine the contribution from unresolved scales to the effective 
flame propagation velocity. However, as a lower cut-off they 
use the laminar burning velocity si. On the basis of the per- 
formed simulations we suggest a slightly higher value in order 
to take into account the cellular structure of the flame below 
the Gibson scale. This may have impact on the nucleosynthesis 
yields from SN la models and will be addressed in a forthcom- 
ing study. 
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